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The spatial structure of single photons pUS] is 
becoming an extensively explored resource used 
for facilitating the free-space quantum key dis¬ 
tribution [4H3 quantum computation [8] as 

well as for benchmarking the limits of quantum 
entanglement generation [3] with orbital angu¬ 
lar momentum modes mm or reduction of the 
photon free-space propagation speed jlQ| . Albeit 
nowadays an accurate tailoring of photon’s spa¬ 
tial structure is routinely performed using meth¬ 
ods employed for shaping classical optical beams 
m uni E], the reciprocal problem of retriev¬ 
ing the spatial phase-amplitude structure of an 
unknown single photon cannot be solved using 
complimentary classical holography techniques 
na HD exhibiting excellent interferometric pre¬ 
cision. Here we introduce a method to record a 
hologram of a single photon (HSP) probed by an¬ 
other reference photon, based on essentially dif¬ 
ferent concept of quantum interference between 
two-photon probability amplitudes. Similarly to 
classical holograms, HSP encodes full information 
about photon’s “shape”, i.e. its quantum wave- 
function whose local amplitude and phase are re¬ 
trieved in the demonstrated experiment. 

The complete characterization of a quantum wavefunc- 
tion of an unknown photon presents a challenging task 
whose difficulty lies particularly in the retrieval of its 
local phase variations. This is caused by the fundamen¬ 
tal property of single photons i. e. their entirely inde¬ 
terminate global phase following from the perfect rota¬ 
tional symmetry of their Wigner functions in the phase 
space [14], which precludes the application of interfer¬ 
ometric techniques such as optical holography utilizing 
fixed phase relation between investigated and reference 
light. Therefore the characterization of photon’s spatial 
structure has never benefited from the precision and a 
simplicity provided by the holography methods nans] 
but as yet has been tackled only using indirect tomo¬ 
graphic techniques [H] or weak values measurements |2]. 

In this paper, we experimentally show that the holo¬ 
gram of a single photon (HSP) encoding full informa¬ 
tion about its spatial structure given by the quantum 
wavefunction 'ip{x) = {x\fj) can be recorded if the 
first-order interference of optical fields is replaced by the 
non-classical interference of spatially varying two-photon 
probability amplitudes. The idea of HSP, sketched in 
Fig. la, relies on overlapping the unknown photon \2pu) of 
an arbitrary local phase profile (p{x) = 8iTg{'ipu{x)) with 


a reference photon \2pr) having the constant local phase 
profile on a beam splitter, both photons occupying sim¬ 
ilar spectral (temporal) modes. Afterwards we measure 
positions of photons which coincidentally left two distinct 
output ports of the beam splitter parametrized by x and 
x' coordinates. Any feature distinguishing photons, such 
as local difference between their quantum wavefunctions 
'ipu{x) and 'ipr{x) prevents them from ideal two-photon 
coalescence known as Hong-Ou-Mandel effect m, thus 
the observation of spatially localized coincidences {x,x') 
serves as a sensitive probe of the spatial structure of the 
unknown photon. As we visualize in Fig. lb, such a co¬ 
incidence event can originate either from transmission or 
reflection of both photons at the beam splitter. These two 




Figure 1. Quantum interference of two spatially struc¬ 
tured photons, a. In analogy to classical holography we 
repeatedly overlap an unknown photon with a reference 
(known) photon |'0r) with the constant local phase profile on 
a 50/50 beam splitter and we spatially localize coincidence 
events in x and x', measuring their joint probability distribu¬ 
tion |T(x,x')|^ which is sensitive to any differences between 
the quantum wavefunctions of the photons 'ipuix) and 
including the local variations of their phases, b. The spatially 
localized coincidence events [x,x) originate from the non¬ 
destructive interference of probability amplitudes of two clas¬ 
sically exclusive, but quantum-mechanically coexisting sce¬ 
narios: (left) the unknown photon in x and the reference pho¬ 
ton in x' have passed through the beam splitter, (right) both 
photons localized conversely in x' and x have been reflected 
from the beam splitter. 
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Figure 2. Encoding of the local phase of quantum wavefunction in the hologram of a single photon (HSP). The 

HSP emerging from a joint probability distribution of the coincidence events encodes (see equation (1^) the local 

phase profile of the unknown photon (f{x). a-c, To illuminate this feature originating from the local phase sensitivity of the 
quantum interference, we depict the expected HSP structure (false colors denotes computed probability gradations) for photons 
in two gaussian modes with identical amplitudes I'lpuix)] = \^r{x)\, differing by the local phase profile of the unknown photon 
presented in the upper-right corner of each plot, a, For the experimental demonstration the purely quadratic local phase profile 
has been chosen, b, HSP for the generic fourth-order polynomial local phase profile, c, HSP for the non-polynomial local phase 
profile resembling a fragment of the Warsaw skyline. 


fundamentally indistinguishable events account simulta¬ 
neously to a two-photon probability amplitude ^(x,x') 
describing one photon localized at position x and the 
other at x', which can be expressed in Feynman’s path 
integral formalism as: 

= ^{{x\ipu){x'\ipr) - {x\ipr){x'\ipu))- (1) 

Thanks to the recent advances in spatially resolved de¬ 
tection techniques such as fiber-coupled avalanche photo¬ 
diodes iniiiH] or single photon sensitive intensified cam¬ 
eras nano] we were able to measure the joint probability 
distribution |^(x,x')p with the resolution high enough 
to reveal its spatial variations originating from the non¬ 
destructive interference of the unknown and the reference 
photon quantum paths. Remarkably this joint probability 
distribution provides information about the local phase 
profile of the unknown photon (p{x): 


\'^{x,x')\^ = i(|V’„(a:;)P|V’r(a:')P + IV’r(a:;)PlV’«(a^')n- 

'^\i)u{x)\\i}u{x')\\i}r{x)\\i)r{x')\cOs{<p{x) -<p{x')). (2) 

The HSP given by |^(x,x')p is entirely insensitive to 
any constant offset of the local phase profile of the un¬ 
known photon in contrast to optical holograms which, 
while being recorded, are extremely sensitive to a phase 
shift between reference and unknown fields. The visibil¬ 
ity of the HSP fringes V is defined by a spectral (tempo¬ 
ral) mode overlap which can be high and stable for pho¬ 
tons generated by different sources such as two indepen¬ 
dent spontaneous parametric down-conversion (SPDC) 


sources ED, quantum dots [22] or even dissimilar sources 
[23] . In Fig. 2 we visualize the HSP structures for exem¬ 
plary local phase profiles. We can resort to one of the 



Figure 3. Experimental setup for measuring HSP. Or¬ 
thogonally polarized unknown and the reference photon, gen¬ 
erated in spontaneous parametric down conversion process, 
are prepared in the same spectral mode. The photons are 
transmitted through the single mode fiber (SMF), separated 
by the polarization beam splitter (PBS) and then, at the out¬ 
put beam waist, the local phase profile ^{x) is imprinted 
on the unknown photon during its double pass propaga¬ 
tion through a phase mask (a cylindrical lens (CLl) for the 
quadratic phase as in Fig. 2a). We localized photons out¬ 
going from two distinct ports of a beam-splitter, here im¬ 
plemented collinearly as a half wave plate (A/2) and calcite 
crystal, by means of the state-of-the-art intensified sCMOS 
camera [191 120 ] • Both the beam waist surface of the reference 
photon and the phase mask surface were mapped onto the 
camera with a phase-preserving 4/ system consisting of two 
spherical lenses (SL). 
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Figure 4. Measured and reconstructed HSP along with the full retrieval of encoded quantum wavefunction. 

a, Directly measured joint probability distribution form an empirical HSP. b, Utilizing independently measured, 

nearly identical amplitudes of the quantum wavefunctions |'0u(a:)|, |'0r(a:)| we numerically reconstructed HSP best matching 
the raw experimental data which closely resembles the theoretically predicted pattern presented in Fig. 2a. c, The measurements 
followed by the numerical reconstruction yields the complex quantum wavefunction of the unknown photon 'ipu (x ), in particular 
its phase (p{x) = arg('0^i(x)). Uncertainty ranges stand for one standard deviation (see Methods for details). 


numerous methods of phase retrieval [24] to infer (p{x) 
from equation as the detection probability distribu¬ 
tions are directly measurable quanti¬ 

ties. 

We selected for experimental demonstration of HSP 
the situation depicted in Fig. 2a where the unknown pho¬ 
ton has the quadratic local phase ip{x) = kx‘^/2R result¬ 
ing in a cross-shaped HSP. Here R stands for a radius 
of curvature, and k = 27r/800 nm for the wave number. 
Both the unknown and the reference photons were gener¬ 
ated via type-II spontaneous parametric down conversion 
process realized in a periodically poled KTP (PPKTP) 
nonlinear crystal pumped with 400 nm light from a con¬ 
tinuous wave diode laser. We ensured the high indis- 
tinguishability of their spectral (temporal) modes, con¬ 
firmed in an independent HOM dip measurement yielding 
the visibility of 91%. 

As presented in Fig. 3, the photons were spatially fil¬ 
tered by a single mode fiber and then separated by polar¬ 
ization beam splitter and directed separately to two arms 
of a delay line and phase imprinting system. The lengths 
of the arms were adjusted to overlap photons temporally 
and to set the constant-phase waists of the mode coming 
out from a fiber collimator on the mirrors surfaces. We 
inserted a cylindrical lens (/c = 75 mm) in the proximity 
of one of the mirrors, thus imprinting the quadratic local 
phase profile in a horizontal direction on the unknown 
photon during its back and forth propagation. Since the 
reference and unknown photon propagating through dif¬ 
ferent arms of the delay line were orthogonally polarized. 


no interference occured at the delay line output. 

The key part of the setup was an intensified comple¬ 
mentary metal-oxide semiconductor (I-sCMOS) camera 
system of parameters suitable to detect spatially resolved 
photon pairs (see refs HU [20] and Methods for camera 
operation details). We imaged the delay-line mirror sur¬ 
faces on the camera using a 4/ system preserving both 
amplitude and the phase of impinging photons spatial 
wavefunctions. A cylindrical lens (CL2) placed in front of 
the camera reduced the mode size in the vertical direction 
perpendicular to the plane of the setup and consequently 
frame reading time. 

The beam splitter transformation was implemented in 
the collinear configuration as a half wave plate followed 
by calcite polarization displacer such that its two out¬ 
put ports corresponded to the two distinct regions of the 
camera. In the experiment we retained for analysis frames 
containing two detected photons, registering their posi¬ 
tions in a horizontal dimension parametrized by x^x' co¬ 
ordinates in respective regions of sCMOS sensor. High 
spatial resolution allowed us to record the subtle varia¬ 
tions of the detected photons positions and thus directly 
measure the empirical coincidence probability distribu¬ 
tion |^(x,x')p. 

The measured HSP consisting of approximately 2.2 x 
10^ detected photon pairs is presented in Fig. 4a, which 
closely resembles the theoretically predicted cross-like 
shape shown in Fig. 2a. Following equation (© we de¬ 
coded the phase (^{x) using one of a numerical methods 
(see ref m and Methods for details), which finds the lo- 
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cal phase profile that yields the coincidence probability 
distribution closest to the measured data as displayed 
in Fig. 4b. The procedure was fed with virtually iden¬ 
tical wavefunctions amplitudes of the unknown photon 
|'0^i(x)|, presented in Fig. 4c, and the reference photon 
\'ipr{x)\ measured independently using the coincidence 
imaging scheme (see ref m and Methods for details). 

We show the complex quantum wavefunction of an un¬ 
known photon i.e. its measured amplitude and the phase 
extracted from its HSP along with the uncertainty ranges 
in Fig. 4c. We found the radius of curvature of the recon¬ 
structed local phase profile of the photon i? = 34 ± 1.5 
mm to be in a good agreement with the value expected 
from a double pass through the phase imprinting lens 
which has been confirmed in an independent measure¬ 
ment by interfering a classical beams in this setup. The 
uncertainty of the reconstructed phase is below 27r/25 
in a central region and it diverges only on the edges of 
the wavefunction due to the scarcity of registered counts 
outside the central region. 

The HSP method naturally transfers the hologram 
recording techniques into the field of quantum optics 
presenting a compelling and promising way of the quan¬ 
tum wavefunction retrieval. The technique can be read¬ 
ily adapted to the more general configurations where 
a reference photon has unknown structure by spatially 
shearing [25] the photons in the second measurement 
run (see Supplementary Information for details). HSP 
technique can be also extended to the two-dimensional 
case requiring the efficient detection of a four dimen¬ 
sional coincidence probability distribution (see Supple¬ 
mentary Information for details). Parallel development of 
low-jitter, time-resolving detectors would allow to readily 
implement HSP in the mathematically equivalent spec¬ 
tral (temporal) domain where local phase sensitivity of 
the non-classical interference has been observed [26] , and 
several wavefunction reconstruction techniques have been 
presented [27]-[29]. Finally let us emphasize that since our 
scheme relies solely on multiparticle bosonic interference, 
it can be generalized for all bosons. Prospective mea¬ 
surement of the hologram of a single atom and further 
retrieval of its wavefunction could utilize the scheme re¬ 
cently reported in the first experimental realization of 
two-boson interference m relying on a similar detection 
technique. 

Methods 

Photon source. The photon pairs consisting of the unknown 
and the reference photon were generated via type-II degenerate 
spontaneous parametric down-conversion (SPDC) process realized 
in 5-mm long periodically poled KTP crystal (poling period 9.2 /rm) 
pumped with 8 mW of 400 nm light from a single mode, continuous 
wave diode laser. The temperature of the crystal was stabilized to 
24.1°C to ensure maximal and stable overlap between the spectral 
modes of generated photons. The photons were spectrally filtered 
by a narrowband 3-nm full width at half-maximum interference fil¬ 


ter, spatially filtered by a single mode fiber, and temporally over¬ 
lapped after polarization beam splitter by means of an optical delay 
line where double pass through quarter wave plate (A/4) rotated 
photons’ polarization by 90°. 

We characterized the indistinguishability of photons used in the 
experiment with the standard avalanche photodiode coincidence 
system by measuring the Hong-Ou-Mandel dip yielding the visibil¬ 
ity of 91%. 

Single photon localization with I-sCMOS camera. To lo¬ 
calize photons with high spatial resolution we used a camera system 
consisting of scientific complementary metal-oxide semiconductor 
endowed with image intensifier camera (I-sCMOS), assembled in 
our group. The image intensifier begins with gallium arsenide pho¬ 
tocathode converting the impinging photons into electrons with 
the quantum efficiency of 20%. Afterwards each electron enters the 
multichannel plate where it triggers the growing charge avalanche 
which hits the phosphor screen resulting in a bright green-light 
flash of a decay time below 200 ns. A typical phosphor flash has a 
diameter of 66 jim and highly random brightness determined by the 
stochastic avalanche process. The flashes are imaged on the sCMOS 
camera sensor via a bright relay lens and real-time localized by a 
software algorithm, which retrieves central positions of the flashes 
from a raw image with subpixel accuracy. We acquired a data from 
1000 X 20 pixel region of interest selected on the sCMOS camera 
sensor corresponding to approximately 11.5 X 10"^ microchannels 
with a frame rate of 7 kHz. We set the time gate of the image in¬ 
tensifier to 30 ns ensuring that virtually no accidental coincidences 
or more than two photons per frame were detected. Moreover, the 
chosen time gate corresponded to a dark count rate of 4 x 10“^ per 
microchannel which could be neglected in further analysis. See Sup¬ 
plementary Information for further details of camera construction 
and operation. 

Measurement details. We measured the amplitude of the 
wavefunctions \'ipu,r{x)\ by setting the half wave plate A/2 (HWP) 
to 0 = 0° and 0 = 45° to direct photons into different output ports 
of the calcite displacer. The nearly identical squared amplitudes of 
both photons were recovered by directly following our coincidence 
imaging scheme m- Then we proceed to the HSP measurement by 
setting HWP to 6 = 22.5° interchangeably measuring wavefunction 
amplitudes and HSP by rotating HWP after each 5 x 10® frames 
out of 1.8 X 10® of total number of collected frames. 

Phase-retrival algorithm. We found efficient to apply a nu¬ 
merical search of ^{x) that matches the measured HSP best, by 
solving the optimization problem according to the general idea sug¬ 
gested by ED- We performed the optimization procedure: 

min I ||T(x,a:')|^ - (x, 11, 

(p(x) 

where |T(x,a:')p stands for the measured empirical distribution, 
(x, is a functional defined by equation (2) constructed 

from the measured amplitudes \'ipu,r{x)\ depending on the vector¬ 
ized phase profile ^{x) to be found, and || • || is the Frobenius norm 
of the matrix. Since the general global search is a computation¬ 
ally hard problem, we divided our optimization into two simpler 
subsequent steps. We assumed that (p{x) is a general fourth-order 
polynomial and we ran a global search finding its coefficient and 
the visibility parameter V. Afterwards we performed the local opti¬ 
mization with unconstrained values of the discretized ^{x), starting 
from the result obtained using the global search. 

Results uncertainties. To account for the uncertainty of the 
empirical HSP in the phase retrieval procedure, we applied the 
Monte-Carlo approach. We repeated the phase retrieval procedure 
5000 times, each time randomizing ^ \'il)u{x)\ and \'il^r{x)\ 

by drawing the initial counts values at each pixel from the corre¬ 
sponding Poissonian distributions. In each realization we obtained 
\^Yec{x^x')\‘^ and the corresponding vector of phase. The Monte- 
Carlo approach resulted in mean reconstructed HSP presented in 
Fig. 4b and the phase profiles whose mean and standard devia¬ 
tion, after unifing their convexities and constant phase offset, are 
presented in Fig. 4c. 
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